Kir form factors and final state interactions in D + — > K 7r + 7r + decays 

D. R. BoitcQ and R. Escribanci 

Grup de Fisica Tedrica and IFAE, Universitat Autbnoma de Barcelona, 
E-08193 Bellaterra (Barcelona), Spain 

Abstract 

We present a model for the decay D + — > K~ir + n + . The weak interaction part of this reaction is described 
using the effective weak Hamiltonian in the factorisation approach. Hadronic final state interactions are 
taken into account through the Kir scalar and vector form factors fulfilling analyticity, unitarity and chiral 
symmetry constraints. The model has only two free parameters that are fixed from experimental branching 
ratios. We show that the modulus and phase of the S wave thus obtained agree nicely with experiment up 
to 1.55 GeV. We perform Monte Carlo simulations to compare the predicted Dalitz plot with experimental 
analyses. Allowing for a global phase difference between the S and P waves of —65°, the Dalitz plot of 
the D + — ► K ~7r + 7r + decay, the Kir invariant mass spectra and the total branching ratio due to S'-wave 
interactions are well reproduced. 
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I. INTRODUCTION 



In 2002, the analysis of D + — ► K ir + ir + decays performed by the E791 collaboration revealed 
that approximately 50% of these decays proceed throngh a low-mass scalar resonance with isospin 
1/2: the 1Tq(800), also called the k [I]. As a matter of fact, the k was the second elusive scalar 
to be firmly detected in D + decays since the scalar-isoscalar /o(600), or a, had been detected by 
the same collaboration in D + — > 7r + 7r~7r + [2]. More recently, the D + —* K~ir + ir + decay was 
revisited by E791 [3| and two other experiments produced analyses based on larger data samples, 
namely FOCUS [4, 5] and CLEO [fj. The main conclusions of the pioneering E791 work have been 
confirmed in both cases. 

In the past, many analyses of Kir scattering data had already claimed the presence of the k pole 
in the scattering amplitude . [lo|. The most precise and model independent determination 

of its position in the second Riemann sheet was produced in Ref. ljj, following the method 
put forward for the a in Using Roy's equations for Kir scattering [13] and Chiral 



Perturbation Theory 
T K = 557 ± 24 MeV 



(ChPT) 



141 ] Descotes-Genon and Moussallam found m K = 658 db 13 MeV and 



Although the experimental results are sound and the k pole is at present theoretically well 
known, a comprehensive and successful description of the reaction D + — ► K~tt + it + is still not 
available (for a recent review see Ref. [3]). Experimentalists, for the want of a better framework, 
commonly fit their data with the isobar model which consists of a weighted sum of Breit-Wigner- 
like propagators. Often, a complex constant is added to the amplitude in order to account for the 
non-resonant decays. It is known, nevertheless, that the adoption of Breit-Wigner functions to 
describe the effect of scalar resonances is problematic. Some of the deficiencies of this approach 
are discussed in Ref. [16] where Oiler pro pos ed the substitution of these functions in the S wave 
by expressions based on unitarised ChPT [171 ] . This model provides a good description of the data 
but, since the weak part of the decay was not tackled, the relative weight of the amplitudes remain 
arbitrary complex parameters to be determined from the fit. 

Little progress has been achieved in the treatment of weak decays of charmed mesons since the 
seminal papers by Bauer, Stech and Wirbel Q]- This fact stems from the mass of the c-quark 
that lies between the heavy and the light domains, rendering heavy-quark approaches or the use of 
chiral symmetry less trustworthy. A first attempt to describe the decay D + — > K~ir + n + from first 

n 

principles was made by Diakonou and Diakonos in Ref. [201. In their work, the weak amplitude 

nn 

was described within naive factorisation with the weak Hamiltonian of Refs. [18|, [19|] and the final 



2 



state interactions (FSIs) were implemented by means of Breit-Wigner type Ktt form factors. They 
considered the contribution of two resonances, namely the iT*(892) and the -^(1430). In the 
light of the present empirical data it is clear that this model cannot provide a good description 



of the decay. In Ref. 



Refs. 



20l | . the decay is mainly driven by the iT*(892) whereas the analyses of 



show that the decay is largely dominated by Kir pairs in an S"-wave state. On 



average, the total scalar signal amounts to 82% 



2l(. Hence, a more comprehensive model for the 



whole scalar contribution is needed to provide a good description of the data. A first step in this 
direction was taken in Refs. 



22 



23J where the irir scalar signal in B —* iririr decays was considered. 



In this framework, factorisation is assumed for the weak amplitude and the ttix scalar form factor, 



constrained by chiral dynamics and unitarity, provides the description of FSIs 24]. In Refs. 



25 



a similar description was utilised to describe the S wave in B — > irirK and B — » KKK decays 



Using the same method, S'-wave FSIs have also been considered in the decay D + — > Tr + ir~ir + [271 ] . 
More recently, Ktt form factors have been employed in the description of FSIs in B^ — ► K tt^tt 
decays 28]. In the present work, we follow the same general scheme where a factorised weak decay 



amplitude is dressed with FSIs by means of non-perturbative Ktt form factors. 

For the weak vertex, we employ the effective weak Hamiltonian of Refs. Q] within naive 



factorisation. Although the assumption of factorisation is less reliable 
has been successfully applied to D decays in several recent papers [27J, 



for the c-quark mass scale, it 



29, 



30 



one should consider the Wilson coefficients as phenomenological parameters to compensate for the 



deficiencies of factorisation 



35] 



The phenomenological values are close to the calculated ones 
but have larger errors than in applications to B decays. The weak amplitude thus obtained 
receives contributions from colour-allowed and colour-suppressed topologies. In the latter, the 
Ktt form factors appear manifestly and the construction of the final state is straightforward. The 
colour-allowed topology is more involved but, assuming the decay to be mediated by resonances 
as suggested 
form factors 
by the Ktt scalar and vector form factors. 

Both form factors have received attention in recent years and are now well known in the energy 
regime relevant to D + — * K~tt + tt + decays. The scalar component was studied in a framework that 



ay the experimental results, the FSIs in this case can also be written in terms of Ktt 
22. 127I . Therefore, in our description the hadronic FSIs are fully taken into account 



36, 



37|, 



. Analyticity, unitarity, chiral 



incorporates all the known theoretical constraints in Refs. 
symmetry, the large-iV c limit of QCD, and the coupling to K-q and Kr( channels were taken into 
account. The results were subsequently updated and we employ in this work the state-of-the-art 
version given in Ref. 



39J. The vector form factor, in its turn, can be studied in r 



KlXV-r 



3 



decays 



40, 



41 



42 



431 ] , where the kinematical range is very similar to the one considered in this pa- 



per. A prediction for this form factor within Resonance Chiral Theory (RChT) 44j] was presented 
in Ref. [^(J and, after the appearance of the detailed spectrum measured by the Belle collabora- 
tion [45I, a fit was performed in Ref. [41]]. Here we employ a slightly different description which 



fulfils analyticity constraints and that was successfully fitted to the Belle spectrum in Ref. 

Our paper is organised as follows. In Section [U we present our model and discuss previous 
treatments of the same decay found in the literature. The numerical results are worked out in 
Section IIIII Finally, we give a summary and discuss the results in Section UVl Details about the 
construction of the Kir form factors employed in this work are relegated to the Appendix. 



II. THEORETICAL FRAMEWORK 

Our phenomenological description of the weak process D + — > K~ir + ir + is based on the effective 
Hamiltonian 

n eS = ^vMidi^o, + c 2 (p)o 2 ] + h.c. , (i) 



where Gf = 1.16637 x 10 5 GeV 2 is the Fermi decay constant 



2l|, V cs V: d = 1 - A 2 , in the 



Wolfenstein parametrisation [46J] with A = sin #,7 = 0.2257 2l|, Ci j2 (/i) are short distance Wilson 
coefficients computed at the renormalisation scale \x = 0(m c ), and 0\ )2 are the local four-quark 
operators 

01 = [Ci7^(l - 7 5 )sip i 7 At (l - yfiuj] , 

02 = [c*7 M (l - l5)Sj][djJ^{l - J 5 )ui\ , 

with = 1, 2, 3) denoting colour indices. At the quark level, the decay D + — > K~ir + ir + is driven 
by the transition c — > sud, i.e. four different quark flavours are involved. In this case, only the two 
tree operators in Eq. ([2]) have to be taken into account. 

The amplitude for D + — > K~ir + ir + is given by the matrix element {K~ir + n + \7i c R \D + ). We 
assume the factorisation approach to hold at leading order (in Aqcd/^-c and a s ) and as a conse- 
quence the amplitude is written in terms of colour allowed and suppressed contributions, Ai and 
A2 respectively, as 

A(D+ -> K-it + it + ) = — £ cos 2 c {aiAi + a 2 A 2 ) + (vr^ vr^) 

v2 

= ^cos 2 0c[ai(^vr+| S V(l - 75)c| J D + )(vr+|^(l - 75 )d|0) (3) 
+a 2 (A-vr+|s7^(l - 75 )d|0)(7r+|n7 /i (l - ^)c\D+)\ + (vr+ ^ vr+) , 



4 



where the last term accounts for the presence of two identical pions in the final state. The QCD 
factors 0-0 are related to Ci^ip) as follows: 



a\(n) = C\(ji) + —C 2 {n) , a 2 {n) = C 2 (//) + — Ci(/x) , 



N r 



Nr. 



where N c = 3 is the number of colours. For these factors we use the phenomenological values 



oi = 1.2 ±0.1 , 



a 2 



-0.5 ±0.1 



(4) 



(5) 



obtained from different analyses of two-body D meson decays 34 ] . 

The non-perturbative hadronic matrix elements in Eq. © involve several Lorentz invariant 
form factors. We first consider those related to the Ai contribution. The transition D + — > K~tt + 
appearing in A\ is more involved and requires a separate analysis. The matrix element from the 
vacuum to the Kir final state is given by 

1 1 11 

F^(q 2 ) + ^ 



(K-ic$\&fd\Q) 



1 ni,< "ViffV) , (6) 



where q = pk + V-k\ an d F+o (q 2 ) are the Kir vector and scalar form factors. Analogously, the 
transition D + — > ir + is given by 



{ir£\u'fc\D'* 



(PD + Pn 2 



2 2 
v/i _ m D ~ m 7r u, 



F?*(q 2 ) ± 



m D — m 



(7) 



where now q = Pd — Pn 2 anci FP^^q 2 ) are the Dtt vector and scalar transition form factors, 
respectively. The amplitude Ai then reads 



An 



m 



Ktt2 



m 



(m\ - m\)(m 2 D - m\) 



7T17T2 
2\^2 



(m K - m^)(m D - m w ) 



mi 



F K *(m 2 K7Tl )F ^(m 2 K 



F¥'"[m 



2 



)F D, { 



m 



Km/ 



(8) 



Dtt/ 2 



ru 



K-K! 



where the Mandelstam variables are defined as 



m Km = (PK + Pm f , rn 2 K = (p K + p„ 



m 2 = fan +P-k 2 ) 2 



with m 2 Kni + m 2 Kn2 + ml 17T2 = m 2 D + m 2 K + 2m\. 

In our analysis, we use a simple pole prescription for the Dtt transition form factors, 



1 — q z m 



pole 



(9) 



(10) 



with m P oie = m^o for the vector case and m po i e = m D *o for the scalar one. The normalisation 
constant is by construction the same in both cases F® n (0) = Fq )7T (0). This parametrisation 



agrees with experiment. The analysis performed by the Belle Coll. on D° — > ir l + v data gives 
for the simple pole model m po i e (l ) = 1.97 ± 0.09 43], which is compatible with the PDG 



value m D *± = 2.01 GeV [2l|. Then, in Eq. JTD]) we take F± n (0) = 0.624 from Ref. [47] and 
m D ,o = 2007 MeV and m D »o = 2.352 ± 0.050 GeV from Ref. [2l|. 

For the Kir vector and scalar form factors, we employ the same expressions that were used in 



431 ] . Since the kinematical region 



the successful reanalysis of r~ — > Kirv T decays performed in Ref. 
for the Kir system available in D — > Kirir decays, txik + m-n < m Kn < mD — m-K, is very similar 
to that of t~ — > Kirv T decays, mx + nhr < m K-n < m r, we consider this choice appropriate. Both 
form factors are constructed such that they fulfil constraints posed by analyticity and unitarity. 
Because of these properties, the form factors satisfy an n-subtracted dispersion relation, which 
in the elastic region admit the well-known Omnes solution 48]. For the Kir vector form factor 



F+ n (s), a good description of the experimental measurement of t~ — > Kirv T was achieved by 
incorporating two vector resonances and working with a three-times-subtracted dispersion relation 
in order to suppress higher-energy contributions [43]. The additionally required scalar Kir form 
factor F^ n (s) had been calculated in the framework of RChT and solving dispersion relations for 



a three-body coupled-channel problem in Ref. 



361 ] . Here, we use the recent numerical update of 



Ref. [39[] . The details of the form factors used in this work can be found in Appendix S 

Now, we turn our attention to the form factors associated with the .Ai contribution. The form 
factor denoting the transition from the vacuum to a pion final state is nothing else than 

(vr^|u7^(l - j 5 )d\0) = ifp W2 , (11) 

where the constant / equals at lowest order in the chiral expansion the pion decay constant / = 
f w = = 130.5 MeV. The form factors related to the transition D + — > K~ir + are more 

complicated. On general grounds, the matrix element (^T~ 7r^" | s'y^ ( 1 — 75)c|D + ) can be written in 
terms of four different form factors j^. But, when saturated with (7r^|u7 At (l — ^y^)d\0) only one 
of those form factors survives, F4, and the amplitude Ai becomes 

Ay = -iUrnlF A (m^ v m 2 K ^) . (12) 

Since this amplitude is proportional to m 2 one would expect it is negligible, as presumed in Ref. 291 ] . 



If this were the case, however, the decay D + — > K ir + ir + would be dominated by the P-wave 



contribution (as demonstrated in Table [nil of Section in contradiction with experiment 2 ll ]. 
This fact forces one to consider the A\ contribution in detail. Unfortunately, the contribution of 
F4 to semileptonic decays, D + — » K~tt + I + vi {I = e,^), is proportional to the lepton masses and 
neglected [5(J |. Consequently, one has to resort to theoretical models. 



Several methods have been considered in the literature. Most of them are based on the as- 
sumption that the D + — > K~ir + transition is driven by intermediate resonances, mainly vectors 
and scalars in this case. We will not take into account the contribution of tensor resonances. In 
the simplest case, one can consider the exchange of a single vector and scalar resonance using 
a Breit-Wigner parametrisation. For instance, in the paper by Diakonou and Diakonos [20] the 
colour allowed contribution is written via the exchange of K*(892) and -fCg(1430) resonances as 



Ax 



^ {K-ir+\K*)(K*\sr(l- 75 )c\D+) | (K-irt\KZ)(KZ\sr(l- l5 )c\D- 



pol "K* 
X ( 7r 2 f |^(l -75)d|0) 



m 



K 



"i 



rn 



A' 



m 







Km 



(13) 



while the colour suppressed contribution is given by Eq. ([8]) but with monopole Kir form factors, 



F*o(q ) = F£J (0)/(l - q /m* ole ), with m po i e = m^ *^ f or the vector and m poh 



m 



for the scalar. Taking the matrix elements from Refs. 



19| one gets 



in 



ml 



m 



An 



m 



Kit, 



M$ r A 



A* (1430) 



(14) 



where N(q 2 ) = m 2 D + m 2 K + 2m 2 - 2m 2 11T2 - q 2 - M(q 2 ), M(q 2 ) = (m 2 K - m 2 )(m 2 D - m%)/q 2 , 
9k*Kk(9k*Kit) are dimensionless couplings associated to {K~k^\K*(Kq)) , and F+ K (m 2 ) and 



Fq K ° (m 2 ) are pertinent vector and scalar transition form factors evaluated at q 2 = m 2 



monopole form is assumed, 



?DK* 



?DK* 



(q 2 



(0) 



-,DK, 



q 2 /m 2 



F DR Z(0) 



pole 



q 2 /m 2 



Again, a 



(15) 



pole 



with m po ie = in both cases 

HQ. 

Experimental data collected in Tables HI and [TT1 indicate that the vector contribution to the total 
signal is largely dominated by the exchange of K* (892). Hence, a Breit-Wigner parametrisation 



20j ]. should be a reasonable approximation 



with a single vector resonance, as considered in Ref. 
to the vector induced signal. This is not the case for the scalar one, where the contribution of 
-?Tq(1430) is marginal. Besides, the possible ETq(800) or k and non-resonant contributions are not 
accounted for in Eq. (I14p . Therefore, a more elaborated prescription taking into account the whole 
scalar contribution is mandatory. Here, we follow Ref. 22J and write the colour allowed amplitude 
Ax in terms of the scalar and vector Kir form factors. We briefly summarise the method applied 
to our case. The D + — > K~ir + matrix element is written as 



(K-tt + \s^(1- 15 )c\D + ) = {K-7T+\R)P R (R\s^(l- 75 )c\D^ 



(16) 



R=S,V 



where we assumed that only scalar and vector intermediate resonances propagate. Tensor reso- 
nances are not included in the sum since the X|(1430) is seen to contribute less than 1% [2lf] . In 
Eq. (fT6l) . (K~ir + \R) is the coupling of Kir to the resonance and Pr stands for the propagation of 
that resonance. The same decomposition is possible for the matrix element which define the scalar 
and vector form factors. Our aim is to substitute the products {K~tt + \R)Pji, usually involving 
Breit-Wigner parametrisations, by expressions based on the relevant form factors. 

For the scalar case, let us take for instance the contribution of 1Tq(1430) alone and write 



-F,fV) = (K-7r + \K* )P R J q 2 )(KZ\sd\0) , 



(K-ir + \sd\0) 

m s -m d 

where the matrix element {Kq\scI\0) defines the scalar decay constant. Then 

1 



(17) 



n^(, 2 ) = (K-«+\KZ)p Rs tf) = o) ms _ md 



(18) 



with XR* a P ure number understood as a normalisation. Hence, the contribution to the matrix 



element in Eq. (I16H is 



<if-7r + |^(l- 75 )c|£> + )| > = X ~F *V)<tfoW(l-75)c|£ + ) . 



\K, 



(19) 



f. Q 



In order to make contact with Ref. [20J one can consider the function Ii R * KlT in a Breit-Wigner 
parametrisation, 



9K*Kn m K* 



(20) 



recovering the scalar contribution in Eq. (|14p . For the remaining matrix element we use 



(K* \s^(l- 75 )c\D^ 



(PD+PR 5 r -^v 



m n ~ m K* n^* 1 



+- 



(21) 



with q = Pd — PR*- Finally, we get the _K"o(1430) contribution to A\, 



(K-7r+|s7^(l- 75 ) C |D+) _ <7r+|u 7/1 (l-7 5 )d|0> 

K* 

= UXR S (m 2 D - myF DKS (rnl)F ^(m Kni : 



(22) 



From Eqs. (|18p and (|20p . one can get an estimate of the absolute value of XR* 



XR* 



Fq 1t {w? 



9k^Ktt 



^(4)1^(4)1 



(4.4 ± 2.8) GeV" 1 



(23) 



where the error includes only the uncertainty in gx*K-K an d ^g*- For the numerical values we have 
used g R * Kn = 3.4±1.9, obtained from B(K% -> Kir) = (93±10)%, T R * = 270±80 MeV [21], and 



\ F o*( m2 K*)\ = 2 - 89 from Ref - |3S1- 

If more than one scalar resonance is exchanged then 



(K-7rf\s^(l - l5 )c\D + )\ s {^\u^{l - 75 )d|0) 



=aEW-WWKK) • l21i 

s 

In Eqs. (|22p and (|24p . the scalar resonances are taken on-shell since it is assumed we are in 
the vicinity of these resonances and hence only small energy regions around the resonance poles 
are considered. However, we want to describe the whole Kir invariant mass range. For such a 
description, we propose the following ansatz for the scalar contribution to A\, 

Af = (K-ir+\s^(l- k)c\D+) | s (7r+|lZ7 M (l-75)d|0) 

= [xsF DS (mD] (ml - m* Kn )F*«{m* Kn ) (25) 

5 

where is a new normalisation constant that contains all the form factors and normalisations 
for the scalar resonances. An estimate for Xs^ * s given by 

xf > Xk S Fo R ° (ml) = (5.5 ± 3.5) GeV" 1 , (26) 

where the value F (m^) = 1.24 ± 0.07 is taken from Ref. y. This value, obtained assuming 
that the form factor is saturated by the Z?j~ pole, is consistent with 1.20 ± 0.07 extracted directly 
32]. Since the estimate in Eq. f|26|) is a lower bound, we prefer to leave Xs^ as a 



from D+ -> K^tt 



free parameter of our analysis to be determined from the reported value of B(D + — ► K~tt + tt + ) 21|. 

For the vector case, let us discuss in some detail the contribution of if* (892). On one side, one 
takes the vector current matrix element in Eq. (|6|) and writes 

(K-ir+ls^dlO) = J2(K-Tr + \K*)P R *(q 2 )(K*\s^d\0) 

K ' poi. (27) 

where (K~ir + \K*) = g R , Klv e{q) ■ (p K - Pn), {K*\s^d\0) = -m R ,f k ,e^*(q), with q = p K + p w , 
and X^poi e ^( ( l) e t(l) = ~9fj.y + Q^v/fn 2 ^. We have made explicit only the contribution of the 
vector transverse degrees of freedom. The dots stand for the longitudinal degrees of freedom which 
can be shown to contribute to both the scalar and vector form factors. However, for the sake of 



9 



comparison, it is enough to consider the transverse part. Comparing Eqs. © and ([27]) . one finds 
the equality 

V-k*kM 2 ) = 9k*k^k*Pk^ 2 ) = - *R*F? n (q 2 ) • (28) 

Ik* 

The former equality must be understood as a replacement of the if* (892) contribution by the 
vector form factor. This replacement should be valid at least in the region around the resonance. 
A direct estimate of Xk* = (4-9 ± 0.2) GeV -1 is obtained using f R * = (205 ± 6) MeV from 



B(t~ -> K*(892)~v T ) = (1.20 ± 0.07)% [2l(]. On the other side, one has 



q»(K*\s^(l ~ ls)c\D + ) = i(e* ■ q)2m R *F»" (^) , (29) 

where the matrix element {K*\s^f^{\ — r y§)c\D + ) is written in general in terms of four different 
form factors [18]. However, after contraction with q = p£> — p R * only the scalar form factor F®^* 
remains 1 . Finally, the if* (892) contribution to A\ is written as 

(K-irf\s^(l - l5 )c\D+)\_^(ir+\u^(l — T5 )d|0> 

= ^(^-^|K*)P^(mLj(^1^(l-7 5 )c| J D+)^ 2 + |n 7M (l-75)d|0) (30 ) 

pol. 

= Un RtK ^mljN(m\,)F^{ml) = UxK*N{m%)F° R * K)F^(m^) . 
In the Breit-Wigner parametrisation, the function Tl R * KlT corresponds to 

n|^U = a m f 9R ^ , (3i) 

m Rsf - m z K7ri - im R *T R * 



again recovering the vector contribution in Eq. (|14p . 

Considering the exchange of more than one vector resonance Eq. (|30f) turns into 



{K-tt+\s^(1 - ^)c\D + ) \ v {tt+\u 1ix (1 - 7s)d|0) 



(32) 



r 

Analogously to the scalar case, we propose to take for the vector contribution to Ax, 
A\ = (^vr+|^(l- 75 ) C |Z)+)| i/ (7r 2 + |n7, t (l-75)d|0) 

= U E [xy^ y K)] iV^J^^kJ (33) 

= UxfN(m Kni )F^(m 2 Kwi ) . 



1 In the notation of Ref. [Til, F + K * corresponds to A$ K 

10 



A lower bound for x'v ls obtained as 

xf > XR*F? R *«) = (4-6 ± 0.9) GeV" 1 , (34) 

where the error takes into account the different results for F® K * (m 2 ) — F® K *(0 ) ex tracted from 
recent analyses. The value F+ K (0) = 0.76 is found in a quark model calculation [511 ] and a lattice 



simulation 52]. This value contrasts with FP K (0) = 1.12 found in Ref. [53J using limits of large 
energy effective theory and heavy quark effective theory. In any case, we like better to leave x'v as 
a second free parameter to be fixed from the experimental value of B{D + — * Kq(892)tt + )+B(D + — > 
i? *(1680)^+) [2l|. 

In Section. 1111} we perform a rather exhaustive numerical analysis of our model and the models 
. For the sake of clarity, our model is defined by the amplitude in Eq. ([3]) resulting 



of Refs. 



20 



from the sum of colour-allowed scalar and vector contributions, Eqs. (|25j) and (|33j) respectively, 

Ai = Uxf(m 2 D - m 2 Km )F K *(m 2 Km ) + f nX f N(m 2 Kni )F^(m 2 Kni ) , (35) 

and the colour suppressed contribution A2 in Eq. ([8]). It is worth mentioning, however, that the 
total scalar amplitude of our model must be re-phased by some amount in order to carry out a fair 
comparison with experimental results, see Eq. (|37p for details. We denote this model as our final 
model. 

III. NUMERICAL RESULTS 

In this section we shall collect all the numerical results arising from the models discussed in 
the previous section and compare them with the experimental results available. Concerning the 
branching ratios, we shall take the PDG averages [2^ shown as the second column of Table HI From 
this table we learn that (i) the contribution of the K~tt + it + mode to D + decays is important and 
accounts for about 10% of these decays, (ii) the decay is strongly dominated by (K~ir + ) pairs in 
the S wave, (iii) although less important the vector K* (892) also gives a sizable contribution, and 
(iv) the branching ratios of submodes containing the next vector and the tensor resonances are 
fairly small. 

Often, the branching ratios for the submodes are estimated from the experimental fit to the 
Dalitz plot through fit fractions. These fractions quantify the weight of the i-th component of the 
amplitude to the final result as 

, = Iv dm K^ dm Kn 2 \Ai\ 2 

Iv dm KK 1 dm2 K* 2 \Y. j A j \ 2 ' 

11 



TABLE I: World average for the relevant branching ratios as reported by the Particle Data Group [2l| and 
results for the three models discussed in SectionflT] {K~n + )s,p denote K~ir + pairs in S or P wave. 
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iviou.61 irorn ivei. 


.4-2 only 
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O.l / /c 


Fixed 


/7 — > 




^ f .04 zt U.ZD J /C 




n 97 Q£ 

U.Z ( /o 


y( . zt \j.z ) /c 


/7 — ? 


l^xv 7T ) P 7T 






z.o4 /o 


Fixed 


D+ -» 


i^ *(1430)vr+ 




0.016% 






,D+ -> 


if*(892) 7r+ 


(1.22 ± 0.09) % 


0.5% 






,0+ -» 


Jf*(1680)vr+ 


(0.16 ± 0.06) % 








,0+ -» 


if 2 *(1430) 7r+ 


(0.030 ± 0.008) % 









In this formula i represents a submode that can be a resonance or the sum of an entire partial wave 



and T> denotes that the integrals are to be evaluated over the whole Dalitz plot (see Ref. 2l|). The 



fit fractions from the analyses of Refs 



are shown in Table |TT] along with the results 



of our model, discussed in the remainder of this section. Experimental groups have used different 
models to fit the Dalitz plot. In Ref. [l| the isobar model was employed and the contribution from 
the k was included as a Breit-Wigner function. In Ref. [4] a isT-matrix model was used for the S 

nriQ 

wave. The results from Refs. [3|, y, [g] are obtained using a quasi-model-independent bin-by-bin 
analysis for the S wave introduced in Ref. Q]. In Ref. Q], a (7r + 7r + )/ = 2 amplitude is also included 
in the model and is found to give a sizable contribution. 

Finally, a comprehensive account of the decay should be able to reproduce not only the known 
branching ratios and fit fractions but also the detailed shape of the Dalitz plot. This is discussed 
for our final model at the end of this section. 



Previous models in the literature 



Here, we update the results of two models found in the literature for the decay under study 



20, 



29] . In both cases the description of the weak decay is based on the effective weak Hamiltonian and 
therefore it is simple to make contact with our model. We begin considering the model presented 
in Ref. [20]] where A± is described by Eq. ([Til) and A2 is given by Eq. (JSj) - We have updated 
the values of the relevant constants for the form factors and for the Breit-Wigner parameters as 
compared with the original work and calculated branching ratios and fit fractions from this model. 
The outcome of this exercise is shown in Tables U and IIIIl The total branching ratio obtained is 
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TABLE II: Fit fractions (in %) for the different submodes of the decay D + — > K~tt + tt + . From Ref. [6fl we 
quote the values for the quasi-model-independent analysis given in their Table VII. Results marked with an 
asterisk are the sum of all contributions to a given partial wave. They do not take into account interference 
effects and were not quoted in the original works. The errors in the results of our model take into account 
the uncertainties in a\ and 02, Eq. ([¥]). 

E791 ('02) [1] E791 ('06) [3] FOCUS ('07) [4] FOCUS ('09) [5] CLEO [6J Our Model 



NR 13 ±7 

k 48 ± 13 

iv *(1430) 12.5 ±1.5 ••• ••• ••• 13.3 ±0.6 

iT*(892) 12.3 ±1.3 11.9 ±2.0 13.6 ±1.0 12.4 ±0.5 9.8 ± 0.5 

iT*(1410) ••• ••• 0.48 ±0.27 

if*(1680) 2.5 ±0.8 1.2 ±1.3 1.9 ±0.8 1.8 ± 0.8 0.20 ±0.12 

#1(1430) 0.5 ±0.2 0.2 ±0.1 0.39 ±0.10 0.58 ±0.12 0.20 ±0.04 

(K-n+) s (73 ±15)* 78.6 ± 2.3 83.2 ± 1.5 80.2 ± 1.4 83.8 ± 3.8 82.0 ±0.3 

(K-7r+)p (14.8 ±1.5)* (13.1 ±2.4)* (16.0 ±1.3)* (14.2 ±0.9)* (10.0 ±0.5)* 15.0 ±0.2 

(7r+7r+) /= 2 ••• ••• ••• ••• 15.5 ±2.8 

J2ifi 88.6 91.9 99.57 94.93 122.8 97.0 



about a factor of 15 smaller than the world average. Moreover, in Table IIIII we show that this 
result is largely dominated by the K*(892) with a fraction of 86.1% of the total result. The 5-wave 
component is represented by the -ffg(1430) alone and accounts for 10.7% of the result. Therefore, 
the model fails to reproduce the absolute branching fractions of Table [J and the strong dominance 
of the S wave that is evident from Tables [J and US As a last comment, note that we employed 
the central values for a\ and a2 given in Eq. (J3J). Shifts within uncertainties in these values could 
produce sizable changes in the branching fractions. However, since the general picture of this model 
does not agree with the known S-wave dominance, we do not attempt to fine-tune these values in 
the case at hand. 

n 

Before turning to our final model it is worth investigating the suggestion of Ref. [29(. In this 
work, the authors advocated that A2 should give the dominant contribution since the colour- 
allowed topology appears multiplied by a factor of ^m^. Following this suggestion, we ignore for 
the moment the colour-allowed topology. In the amplitude Ap. the Ktt form factors enter manifestly 
and it is straightforward to introduce the ones from Refs. 3{], [3], as shown in Eq. ([5]). The use 
of these form factors improves the description of FSIs as compared with Ref. [2^] incorporating 
constraints from analyticity and unitarity. The numerical results from this model are shown TablesU 
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and IIII1 The total branching ratio is now about a factor of 3 smaller than the experimental average. 
Nevertheless, the dominant contribution is again given by the P wave that accounts for 89.6% of 
the total result. Since the result for A2 is unambiguous and we employed state-of-the-art form 
factors we are led to the conclusion that Ai must be taken into account. As a matter of fact, it 
is now transparent that the large 5-wave contribution originates precisely in the colour-allowed 
topology. 

TABLE III: Results for the fit fractions (in %) arising from the three models described in the text. 



Mode Model from Ref. [20] A 2 only Our Model 
NR 

K$(U30) 10.7 

K*{892) 86.1 

{K-ir+) s ■■■ 8.5 82.0 ±0.3 

{K-tt+)p ■■■ 89.6 15.0 ±0.2 

J2ifi 96.8 98.1 97.0 



Our model 



determined in a bin- by-bin basis directly from data. This is done assuming a reference amp 
customarily that of the K*(892). The analysis was repeated by CLEO [6j and FOCUS 



Let us now investigate in detail the numerical results for our final model which includes the 
contribution of both A\ and A2 topologies. The corresponding expressions are given in Eqs. (JSJ) 
and (|35p . We begin by considering the ,5- wave description which is, in our opinion, the main 
aspect of the problem. On the experimental side, in 2006, E791 introduced a new type of Dalitz 
plot analysis [3J where, instead of modelling the Kir S wave, its absolute value and phase are 

itude, 
5J] with 

similar results. It is important to remark that this framework can only be considered as quasi- 
model-independent (QMI) since the P and D waves are still described by their isobar expressions. 
Nevertheless, since the isobar prescription seems to be more accurate for these latter waves, one can 
expect that the results have little model dependence. Therefore, the results of the QMI analyses 
of Refs. jj, 0, f| are the better source of empirical information about the Kit S-wave amplitude in 
D+ -> K-ir + ir + . 

The QMI measurement of the S- wave phase can be used to test whether Watson's theorem 54] 



holds for the three-body decay in question. The theorem states that, in the elastic domain, the 
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Kir S wave would exhibit the corresponding Kir scattering phase shift. However, this is valid only 
in the absence of genuine three-body effects. Therefore, in D + — » K~ir + ir + , the empirical S- wave 
phase could be distorted as compared with the scattering one due to interactions of the resonant 
Kir pair with the bachelor pion. In our model, the S'-wave FSIs are described by the Kir scalar 
form factor of Ref. 



39| in a quasi two-body approach, i.e., we assume that the Kir pairs in Eq. ([3]) 
form an isolated system and do not interact with the bachelor pion. Moreover, the form factor of 
Ref. [39f] is obtained from dispersion relations that fix its phase to be the scattering one within 



the elastic region 36]. Consequently, our S-wave amplitude has the Kir 1 = 1/2 scattering phase 
up to roughly 1.45 GeV where the Krj channel starts playing a role. We compare in Fig. [T]the 
experimental results from Refs. with the phase of our S wave. The high-statistics results of 

CLEO collaboration have the smallest errors. One observes from Fig. [1] that the QMI phases start 
at negative values ranging from -60° |f| to -145° Q] whereas our phase evolves from 0° up to about 
200° (modulo ir) in the allowed phase space. Since we are dealing with a production experiment, 
a global phase difference is expected as compared with scattering results [4j. Therefore, we allow 
for a global phase shift a in our S-wave amplitude 2 : 

As (m 2 K7Ti , m 2 K7T2 ) -> e ia As (m 2 KlTi , m 2 Kn2 ) . (37) 

In Fig. [TJ we also plot as the dot-dashed line the phase of our amplitude shifted by a = —65°. 
With this shift, we see that up to 1.5 GeV CLEO's results and ours share a remarkably similar 
dependence on energy 3 . The results of E791 and FOCUS seem to have a somewhat different 
energy dependence, although they have larger error bars due to smaller statistics. Inspired by the 
inspection of Fig. [TJ we consider as our final model the one given by Eqs. ([8]) and ()35[) with a shift 
of a = —65° in the S-wave phase as defined in Eq. ([371) . We will discuss further consequences of 
this shift below. 

In order to compare the absolute value of our S wave amplitude with experimental data, we 
need fix the only two free parameters that occur in our model, namely the normalisation constants 
and Xv ■ Estimates for the normalisations were given in Eqs. (I26p and (|34|) but in order to 
perform a careful comparison with experimental results we choose to refine these values. With 
that aim, we employ the following strategy. The constant X]/ * s fi xe d in order to reproduce the 



2 The results of the model are sensitive only to the phase difference between the S and P waves. Therefore, a can 
be considered as a global phase difference between the two waves. 

3 The 5- wave phase of the form factor of Ref. [3!j exhibits around 1.8 GeV a deep similar to the one observed in 



the experimental results of Fig. [TJ 



15 



to 
CD 
CD 

1 

O) 

CD 

33 

CD 
in 
cc 



CD 
> 
CC 

5 
CO 



250 



200 



.-, 150 



E791 [3 
Focus [5 
CLEO [6; 
Our phase 
Our phase -65° 



100 




-100 



-150 



-200 



FIG. 1: (colour online). S-wave phases from the QMI analyses of Rcfs. . The solid line is the phase 

of our S'-wave amplitude with a = 0° in Eq. (|37p . whereas the dot-dashed line is the S'-wave phase with 
a = —65°. The dashed line delimits the Krf threshold. 



value of the sum of all vector submodes 4 in the second column of Table HI Then, we fix the scalar 
normalisation x*s requiring the total branching ratio from our model to match the world average 
of Table HI Taking the central values for cq and 02 given in Eq. (JH) this procedure gives 



Xs 



f = 4.9 ± 0.4 GeV -1 



Xv 



4.4 ±0.6 GeV -1 



(38) 



in good agreement with our estimates in Eqs. (|26|) and (|34|) . The uncertainties take into account 
the error in oq and 02 which dominate by far as compared to the relatively small errors of the world 
averages of Table HI We are now in a position to compute the scalar branching ratio, shown in 
Table HI as well as the fit fractions of the total vector and scalar contributions which are shown in 
Tables HI and IIIII The model reproduces the dominant S- wave contribution and gives fit fractions 
in fair agreement with the experimental results. 

We can now compare the absolute value of our 5-wave amplitude with experimental results from 
the QMI analyses. However, since in isobar-like analyses the fit is sensitive only to the relative 
weights of the amplitudes, in order to compare the measurements with our result we need perform 



4 This procedure does not take into account possible interference effects. However, these effects are likely to be small 
since the resonances are relatively narrow. Furthermore, to the best of our knowledge, there is no experimental 
value for the total P-wave branching ratio. 
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a normalisation. We define a normalised 5-wave amplitude by 

As 



Af-{ml^ml n2 ) 



1/2 • 



(39) 



(iv dm K^ dm K^ \M 

This amplitude, by construction, is free of any global constants that appear in *4c and has dimension 



of [Energy] 2 . Interpolating the results from the tables found in Refs. 



flfl 



we can calculate 



the normalised S wave for each experiment. We repeated the same procedure for our total S- wave 
amplitude. The QMI results for the S wave are compared with our model in Fig.[2j Up to 1.55 GeV 
the agreement of our results with the experimental ones is remarkable. 




1.2 
m K3t [GeV] 



FIG. 2: (colour online). Absolute value of the S wave measured in Refs. 0,0,3 compared with our model. 
The amplitudes are normalised according to Eq. (|39|) . The dashed line delimits the Krf threshold. 



Finally, we can perform a Monte Carlo (MC) simulation to obtain a Dalitz plot from our model 
and compare the diagram and its projections with experimental results. For the lack of a true data 
set, we resort to a MC simulation of the original E791 data [l|. Reproducing their fit function, 
we generated a symmetrised Dalitz plot with 14185 independent signal events which corresponds 
to 6% of background contamination in the total sample [1]. The obtained diagram is shown in 
Fig. [3^1. Then we performed the same exercise for our model and the result is shown in Fig. [3b. 
It is important to remark that the shape of the Dalitz plot is related to the global phase shift of 
Eq. (|37p . In the words of Ref. sj], the asymmetry in the Kit P-wave bands reflects the value of 
a. We have checked that taking a = 0° in Eq. (|37|) reverses the observed asymmetry, i.e., the 



17 



high-energy part of the Dalitz is more populated than the low-energy corner. Consequently, we 
confirm the finding of Ref. Q]: the asymmetry pattern in the Dalitz plot is a direct consequence 
of the global phase difference between the S- and P-wave phases. Finally, in Fig. 0] we show the 
projections of the diagrams of Figs. E^, and [3b. The results for our model with a = —65° and the 
simulated E791 data agree quite well. The discrepancy in Fig. around 1 GeV 2 is due to the 
interference pattern between the S and P waves. This could be fixed through a fit to real data, 
which would give a refined value for a. One also sees around 2.5 GeV 2 a second discrepancy, seen 
in both Figs. 0^ and 0b, that is a consequence of the disagreement of our S wave with respect to 
the experimental ones for rriK-n > 1-45 GeV, as shown in Fig. [2j Small isospin-breaking effects in 
the P-wave are to be expected as well, since the vector Kir form factor employed here was obtained 
from t~ — > Kit v t decay data [43|] where the charged vector resonances intervene. 
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FIG. 3: (a) Monte Carlo simulation for the Dalitz plot of the E791 original analysis [l| (b) Same for our 
model with a global shift of —65° degrees in the S'-wave phase (see text and Fig. [T|). The number of 
independent events is 14185, which correspond to the estimate of the signal events in Ref. 



As a final comment, we remark that we do not include (7r + 7r + )/ = 2 interactions in our model. 
Within the framework employed here this contribution does not appear. Since the inclusion of 
an ad- hoc 1 = 2 amplitude would downgrade the model, we prefer to consider only the 1=1/2 
FSIs. Additionally, (7r + 7r + )/=2 scattering is entirely non-resonant [21] with a slow variation of the 
corresponding phase shift 55[], indicating that interactions in this channel are weak. Furthermore, 
from an experimental point of view, the need for the 1 = 2 amplitude is not well established and 
requires further confirmation (see Table ITT]) . 
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FIG. 4: (colour online). Projections from the MC generated Dalitz plots of Figs. [3^ and[3jD. The error bars 
and the bands represent solely statistical fluctuations, (a) Total projection, (b) high-energy projection, (c) 
low-energy projection. 



IV. SUMMARY AND DISCUSSION 

We have presented a model aimed at describing the decay D + — > K~ir + ir + . The weak amplitude 
is described within the effective Hamiltonian framework with the hypothesis of factorisation. The 
Kit hadronic FSIs are treated in a quasi two-body approach by means of the well defined scalar and 
vector Kir form factors, thereby imposing analyticity, unitarity and chiral symmetry constraints. 
We used the experimental values for the total and P-wave branching ratios to fix the two free 
parameters in the model. The relative global phase difference between the S and P waves was 
fixed phenomenologically using the experimental results of Ref . [(J . 

The use of the Kit scalar form factor is shown to provide a good description of the S"-wave 
FSIs. Both the modulus and the phase of our S wave compare well with experimental data up 
to rriK-w ^1-5 GeV. It is worth mentioning that the form factor we used has a pole that can be 
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identified with the k. Furthermore, the model is able to reproduce the experimental fit fractions 
and the total 5-wave branching ratio. Finally, the Dalitz plot arising from the model agrees with 
a MC simulated data set. 

The main hypotheses of our model are the factorisation of the weak decay amplitude and the 
quasi two-body nature of the FSIs. Therefore, the success of our description for mxn ^ 1-5 GeV 
suggests that, in this domain, the physics of the decay is dominated by two-body Kir interactions. 
We are led to conclude that effects not included in our model such as the 1 = 3/2 non-resonant 
Kit S wave, the non-resonant 1 = 2 ir + ir + interactions and genuine three-body interactions, could 
be considered as corrections to the general picture described here. 

Part of the discrepancy observed in our Dalitz plot is due to the disaccord of our S"-wave 
amplitude for mjc-K ^1-5 GeV. A possible cause for this disagreement is the fact that factorisation 
in a three-body decay is expected to break down close to the edges of the Dalitz plot 
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Furthermore, in this region, the kinematical configuration of the final state momenta renders the 
quasi two-body treatment less trustworthy as well. Finally, our model does not include the tensor 
component. Although marginal, this amplitude has a non-trivial distribution in the phase space 
and could induce sizable interference effects in our plots. In the vector channel, we find puzzling 



that the A"* (1410), which gives a sizable contribution for r~ — ► Kirv r 4l|, |43J, is hardly seen in 
experimental analyses of D + — ► K~ir + ir + . 

In conclusion, since we do not fit the Dalitz plot we think that the agreement between the model 
and the experimental data is satisfactory. 
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APPENDIX A: Ktt FORM FACTORS 



The scalar and vector Ktt form factors employed in this work were obtained respectively in 



Ref. [39| and Ref. 43J]. The details can be found in the original references but for the sake of 



completeness we briefly summarise here how they are obtained. 



Scalar Ktt form factor 



The framework for the determination of the scalar Ktt form factor, F^ n (s) : is described in 



detail in Ref. [361 ] . The results were numerically updated later and we employed in our numerical 
analysis the latest version given in Ref. fl. In Ref. Q, the authors solved a generalised Omnes 
problem where three channels, namely Ktt, Kr\ and Krj , are taken into account. In this framework, 
the scalar form factor for channel k, Fq(s) (where 1 = Ktt, 2 = Krj and 3 = Krj'), can be cast as 
a sum over the three channels as 

^.ifL^ww^w, (A1) 

7r f— f J (s — s — ie) 



s, 



In the last equation, Sj is the threshold for channel j, dj(s) are two-body phase-space factors 
and tg~ ^ are partial wave T-matrix elements for the scattering k — > j. The form factors are 
obtained solving the coupled dispersion relations arising from Eq. (IA1|) . This is done imposing chiral 
symmetry constraints and using T-matrix elements from Ref. [10j] that provide a good description 
of scattering data. One recovers the elastic approximation by considering solely the contribution 
of the channel k to the right-hand side of Eq. (jAlj) . which is then reduced to the usual Omnes 
equation 

2. Vector Ktt form factor 



The vector Ktt form factor, F+ n (s), employed in this work was obtained in Ref. 43] within a 
dispersive representation from fits to t~ — > Ktxv t data obtained by the Belle collaboration [3]. 
The reduced vector form factor F+(s) = F^ 7 * (s)/ F+ n (0) is written in terms of a three-times- 
subtracted dispersion relation that takes the form 

Scut 



F?*(s) = exp 



•s 



1 s 2 s 3 f , S^(s') 



a i — 2 + o" 2 — 4 ^ — / ds 'ri 



m 2 2 m% tt J (s') 3 (s' — s — iO) 



(A2) 
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where sk-k is the Ktt threshold and <5f n is the form-factor phase. The subtraction constants ct\ 
of F+(s) around s = 0, 



and «2 can easily be related to the slope parameters A\ , which appear in the Taylor expansion 



F^(s) = 1 + A' + 4 + h'i^ + ■■■ , (A3) 

as A+ = a\ and A" = 02+ af. The cutoff s cu t is introduced as the upper limit of the Omnes integral 
to study the importance of the high-energy region which is strongly suppressed by the factor s' 3 in 
the denominator of the integrand of Eq. (|A2j) . Furthermore, within the elastic region, 5f- n is the 
P-wave 1 = 1/2 Kir scattering phase shift. An advantage of the three-times-subtracted form of 
F+ n (s) is to make the results less sensitive to deficiencies of the phase shift in the higher-energy 
region. Then, the integral in Eq. (|A2j) emphasises the lower-energy domain (elastic domain), for 
which one can provide a reliable model for the phase shift. The description of <5f ^ we used is 
inspired by RChT and includes the contribution of two vector resonances namely the K*(892) and 
the i^*(1410). The detailed expressions can be found in Ref. [431 ] . 
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